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ABSTRACT 

We examine the possible reionization of the intergalactic medium (IGM) by the som-ce UDF033238.7- 
274839.8 (hereafter HUDF-JD2), which was discovered in deep HS T/ YLT / Spitzer images obtained as 
part of the Great Observatory Origins Deep Survey and Hubble Ultra-Deep Field projects. Mobasher 
et al. (2005) have identified HUDF-JD2 as a massive (~ 6 x lO^^M©) post-starburst galaxy at redshift 
z> 6.5. We find that HUDF-JD2 may be capable of reionizing its surrounding region of the Universe, 
starting the process at a redshift as high as z« 15 ± 5. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: intergalactic medium — 

galaxy: individual (HUDF-JD2) — cosmology: observations — cosmology: theory — 
cosmology: early universe 



1. INTRODUCTION 

The reionization of the intergalactic medium (IGM) is 
undoubtedly one of the most crucial events in the evolu- 
tion of the Universe. Since reionization drastically reduced 
the opacity of the Universe to ionizing radiation, it is likely 
to have infiuenced and/or determined the formation and 
evolution of galaxies and other structures in the early Uni- 
verse (see, e.g., Loeb & Barkana 2001). The presence of 
strong Lyman-a forest absorption in the spectra of dis- 
tant quasars indicates that reionization was completed by 
z « 6 (Becker et al. 2001, Fan et al. 2002), while the 
polarization of the cosmic microwave background (CMB) 
radiation suggests it may have begun at a much earlier 
epoch, z ~ 17±5 (Spergel et al. 2003). Theoretical efforts 
to understand and constrain the time evolution of cosmic 
reionization sources {e.g., Loeb & Barkana 2001; Stiavelli, 
Fall & Panagia 2004a, hereafter referred to as SEP) have 
to be confronted with observational identifications of the 
ionization sources. 

The Hubble Space Telescope (HST) observations ob- 
tained by the Great Observatory Origins Deep Survey 
(GOODS) project (Giavalisco et al. 2004), and more re- 
cently the observations of the Hubble Ultra-Deep Field 
(HUDF; Beckwith et al. 2005) have opened a window 
through which one may hope to peek and see distant 
sources that could provide a direct clue to the reioniza- 
tion process. The available data have already extended 
the searches to z ^ 6.5 by the drop-out techniques (Dick- 
inson et al. 2004, Giavalisco et al. 2004, Bunker et al. 2004, 
Yan & Windhorst 2004, Bouwens et al. 2004b, Eyles et 
al. 2005) and narrow-band surveys (Rhoads et al. 2004, 
Stern et al. 2005, Taniguchi et al. 2005). Yet mea- 
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sures of the luminosity functions remain uncertain (Yan & 
Windhorst 2004, Bouwens et al. 2004a), and the identifi- 
cation of the reionization sources is still just a possibility 
[e.g., Stiavelli, Fall & Panagia 2004b). This is because, 
even with the exquisite images of the HUDF, it is hard to 
identify and characterize sources at redshifts appreciably 
higher than 6, which are crucially important to properly 
define and understand the processes associated with reion- 
ization. 

At present, combined ultra-deep images of faint galax- 
ies taken by the Hubble Space Telescope and Spitzer Space 
Telescope (Werner et al. 2004) can be used to provide the 
multi-waveband information needed to search for galaxies 
at even higher redshifts and to explore their nature. In 
particular, the HUDF observations (Beckwith et al. 2005, 
Thompson et al. 2005), representing the deepest images of 
the Universe at optical (B, V, i, and z bands; HST/ACS) 
and near-infrared (J, H; HST/NICMOS) wavelengths, 
combined with the VLT and Spitzer/IRAC observations 
obtained by the GOODS project at NIR (if,; ESO/VLT) 
and mid-infrared (3.6-24 ^m; Spitzer/IRAC-|-MIPS) wave- 
lengths (Vandame et al. 2005, in preparation, Labbe et 
al. 2005, in preparation, Dickinson et al. 2005, in prepa- 
ration), are well suited to accomplish this aim. Analyz- 
ing these data, Mobasher et al. (2005) searched for very 
red galaxies with (J — H)ab > 1-3 and no detection at 
wavelengths shorter than J (i.e. J— band dropouts). They 
found two sources that satisfy the J— band dropout cri- 
teria, UDF033242.8-274809.4 and UDF033238.7-274839.8 
(hereafter HUDF- JDl and HUDF- JD2, respectively). The 
two sources are well detected in the K-band (ISAAC). 
Both sources are clearly detected in all four IRAC channels 
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(3.6, 4.5, 5.8 and 8 fim). 

The first of theses sources, HUDF-JDl, appears to be 
positionally associated (within 1 arcsec) with an X-ray 
source. Since no acceptable fit to its observed spectral 
energy distribution (SED) could be found in terms of stel- 
lar synthesis models, Mobasher et al. (2005) conclude that 
it is most likely to be an obscured AGN at intermediate 
redshift. 

On the other hand, the presence of a clear Balmer break 
in the observed SED of the HUDF-JD2 solidifies its high 
redshift identification and reveals a post-starburst popu- 
lation. From their best-fit models, Mobasher et al. (2005) 
derive a photometric redshift of z ^ 6.5 and a bolometric 
luminosity of Lboi = 1 x 10^^ (for a cosmology with 
Ho = 70 km s-i Mpc^i, = 0.3 and Oa = 0.7). Using 
the M*/Lboi ratio from the model fit, the mass in stars is 
estimated to be M* ~ 6 x 10^^ Mq for an adopted Salpeter 
initial mass function with minimum and maximum stellar 
masses of 0.1 and 100 Mq, respectively. 

The strong Balmer break displayed by HUDF-JD2 re- 
quires that its SED be dominated by stars of spectral types 
AO or later, i. e. main-sequence stars with masses less than 
^ 3Mq {e.g., Allen 1973). This fact unambiguously in- 
dicates a rather old age (>300 Myr) for the stellar pop- 
ulation of HUDF-JD2. Mobasher et al. (2005) conclude 
that the stars were formed at z> 9 {i.e. an age of the Uni- 
verse tu < 540 Myr), and possibly as high as z~12-20, 
so that most of the stars were in place when the Universe 
was only ~ 200 — 400 Myr old. Moreover, their model fit- 
ting indicates that HUDF- JD2 formed the bulk of its stars 
very rapidly, on time scales < 100 Myr, and the subse- 
quent evolution was essentially passive. A stringent upper 
limit to the starburst age is set by the photometric red- 
shift z ~ 6.5 of HUDF-JD2, which corresponds to a time 
when the Universe was only 830 Myr old. Adopting this 
prior, the model fitting procedure indicates that the age of 
the stellar population in HUDF-JD2 is likely to be brack- 
eted between 350 and 650 Myr. These ages correspond to 
redshifts of galaxy formation between 10 and 20. 

The estimated reddening is quite modest, with E{B — 
V) < 0.06 at 95% confidence limit. The best-fit model 
SED also provides estimates of the overall metallicity of 
the galaxy, which is bracketed within the interval Z = 
0.02 - 0.004, i.e. between Zq and Zq/5. In the following, 
we will use a fiducial value of Z = 0.008 = O.AZq, keeping 
in mind that this value may be uncertain by a factor of 
two. 

Given the high redshift of HUDF-JD2 and the short time 

scale over which its stellar population was formed, it is of 
interest to explore whether the radiation from this galaxy 
was able to reionize its surrounding intergalactic medium. 
In this letter, we estimate the output of ionizing radiation 
from HUDF-JD2 over its past history. We will show that 
the Lyman continuum flux radiated by this source in the 
early phases of its evolution may be enough to reionize its 
surrounding IGM appreciably and, under favorable con- 
ditions, even ionize it completely. For consistency with 
Mobasher et al. (2005) paper, we adopt Hq = 70 km s~^ 
Mpc-\ 0„ = 0.3 and Qa = 0.7. 



Reionization is a process that depends on the UV out- 
put of galaxies integrated over the time interval when they 
are active UV sources. In this respect, HUDF-JD2 is a 
promising reionization source in that, according to the 
models, its integrated Lyman continuum photon flux is 
about 4 X lO""^ photons that had to ionize the neutral hy- 
drogen gas within a comoving volume of about 14—18 x 10^ 
Mpc^, the exact volume depending on the specific redshift 
(within the range z ~ 9 — 15; see below) at which H was 
fully ionized. For a comoving hydrogen density of 6.2 x 10*'^ 
atoms Mpc~^ {e.g., Spergel et al. 2003), it appears that 
the produced Lyman continuum photons outnumber the 
H atoms by almost an order of magnitude and that, there- 
fore, a complete reionization is a real possibility. However, 
a realistic treatment of this process should take into ac- 
count the facts that only a fraction / of the Lyman contin- 
uum photons escape from the source, and that the ionized 
gas may rccombinc, and do so at a particularly high rate if 
the IGM material is clumped, traditionally parametrized 
by the clumping factor, C > / < n-H These ef- 

fects tend to reduce the overall efficiency of the ionization 
process. 

Here, we take advantage of the extensive set of model 
calculations made by SEP who considered a wide range of 
semi-analytical models to estimate the reionization of the 
IGM. All cases studied by SEP were based on the implicit 
assumption of a "snapshot" of a statistical assembly of 
cosmic som'ces, which possibly formed at different epochs 
and were efficient UV radiation emitters over different time 
intervals, but, on average, constituted a steady supply of 
ionizing photons over the chosen redshift interval. The 
situation studied here is quite different. HUDF-JD2 is 
a single, powerful source for which both the birth time 
and subsequent evolution are constrained by the observed 
SED. Therefore, we have to allow explicitly for its past 
evolution. 

Although the physical quantity that determines the 
reionization is the total {i.e. time integrated) Lyman con- 
tinuum photon flux, for a convenient comparison with ob- 
servations, SEP parametrize their models in terms of the 
ratio of the Lyman-continuum photon flux to the UV flux 
at 1400A, Cjon = N^yc/ {i'L^)i4oo, which by construction 
remains constant over the time reionization takes place. 
For the case of HUDF-JD2, it is convenient to consider the 
ratio of a time averaged Lyman-continuum photon flux to 
the current UV flux at 1400A. First, we define the average 
Lyman-continuum photon flux as: 

[! NLyc{t)dt 

< NLyc{t) >- . ; ; X ' T 

min[(t - ti),Atreion\ 

where NLyc{t) is the instantaneous Lyman-continuum 
photon flux of the galaxy, ti is the time of the galaxy birth, 
and Atreion is the time interval during which reionization 
takes place. We adopt Atreion = 100 Myr, which coincides 
with the approximate duration of the starburst in HUDF- 
JD2 as estimated by Mobasher et aL(2005). We can relate 
the average Lyman-continuum photon flux < NLyc{t) > 
to the luminosity in the detectable UV by defining an ef- 
fective ionization parameter: 



2. IMPLICATIONS FOR REIONIZATION 
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Figure 1 shows the ionization parameter as a function 
of age for a starburst lasting 100 Myr, which is com- 
puted using Starburst99 models (Leitherer et al. 1999) 
for different metallicities and adopting a Salpeter initial 
mass function within the interval 1 and 100 Mq. For 
our range of metallicity, the value of tion at about 350- 
650 Myrs (i.e. the range within which the HUDF-JD2 age 
is likely to lie) is rather high. For our fiducial metallic- 
ity of 0.4Zq = 0.008 and for an age of 500 Myr, we find 

^lon — ^ ^^^^ erg~^, i.e. 4.8 times higher than the value 
for a 100,000K blackbody source, even though in the syn- 
thesis models no star ever attains such a high effective 
temperature. This is because the ionizing flux was very 
high at the early epochs, while the observed 1400A flux 
correponds to a later epoch when the starburst has weak- 
ened considerably. We note that since the ionization pa- 
rameter depends rather strongly on the age and metallicity 
of the galaxy, e^// may be uncertain by a factor of two. 
The results are nearly the same for any starburst duration 
shorter than 100 Myr. In this context, we note that a rela- 
tively long starburst (i.e. duration '^lOO Myr) is probably 
needed to account for the high metallicities implied by the 
best-fitting SED and that these, in turn, may require mul- 
tiple episodes of nuclear reprocessing. 

Scaling the results of SFP model calculations for 10^ K 
(their Table 1, and Appendix A), we have computed the 
constraints on source surface brightness for a case in which 
reionization starts at = 13 and is complete hy z = 10.7 
(i.e. about 500 and 400 Myr, respectively, before the epoch 
corresponding to our observations) with ionization effi- 
ciency elll = 9.6 X 10^^ erg~^ (see Figure 2). In this 
diagram the scaling factor is simply the ratio of the e^£{ 
value to the reference value of SFP for Tbb = 100, OOOK, 
eio„(100,000ii:) = 2 X lO" erg-^ The one difference rela- 
tive to SFP model calculations is that here the only contri- 
bution to the near UV flux is from stellar radiation because 
the nebular emission of the galaxy is negligible at the time 
it is observed. This simple scaling procedure can be em- 
ployed to describe all cases involving individual sources 
with well defined evolutionary histories. 

In our case, we find that with the observed magnitude 
at A ~ 1.1/xm (corresponding to rest-frame wavelength 
1400A), ruAB — 27.0, and an effective source density of 
N ~ 1/(2.5)^ = 0.16 arcmin"^ (i.e. one source within 
the NICMOS field, 2.5' x 2.5'), HUDF-JD2, as an isolated 
source, could have ionized the hydrogen gas within the vol- 
ume encompassed within the solid angle subtended by the 
HUDF-NICMOS boundaries up to the redshift z ~ 10.7 
(i.e. the "cell of Universe" that we can associate to it) 
only if the escape fraction of Lyman-continuum photons is 
/ > 0.26 (for C=l) or the clumping factor is C < 21 (for 
/ = 1) (see Figure 2). This is already an impressive result 
because we are dealing with one source, which, even with 
rather severe assumptions, such as / ~ 0.2 and C ~ 20, 
can be responsible for at least 20% of the ionization of its 
Universe cell. 

By taking an escape fraction less than unity we are im- 
plicitly including the effects of absorption of Lyman con- 
tinuum photons by both neutral hydrogen and dust grains 
within HUDF-.JD2. However, in our calculations we also 
have to allow for dust extinction that may attenuate the 



source spectrum at wavelengths longer than the Lyman 
hmit. The best-fit to HUDF-JD2 SED yields an estimate 
of a rather modest selective extinction E{B — V) < 0.06 
(Mobasher et al. 2005), which implies an upper limit to 
the extinction at rest frame 1400A of about 0.54 magni- 
tudes. This result is retatively independent of the adopted 
extinction law, be it the effective attenuation law for star- 
burts galaxies {e.g., Calzetti et al. 1996, Chariot & Fall 
2000), or the recent measurements of the extinction law 
in a galaxy at 2: = 0.83 (Mediavilla et al. 2005). The ex- 
tinction correction increases the source flux at rest frame 
1400A by a factor of ~ ioO-4xo.54/2 ^ 28 on average, 
thus moving the source to the right in Figure 2 (open tri- 
angle), indicating that for / ~ 0.2 and C ~ 20HUDF-JD2 
by itself could account for up to 25% of the total ionization 
of its surrounding IGM. 

In addition, IIUDF-JD2 may not be an isolated source, 
but rather the brightest member of a population of galax- 
ies within one and the same volume, whose other members 
are too faint to be detected in the available images. Adopt- 
ing a Schechter luminosity function with a slope a = 1.6 
and a cutoff magnitude M^jy = —21.2 (as appropriate for 
z=3 galaxies; Stcidel et al. 1999), with the condition that 
HUDF-JD2 is the only source brighter than thab — 27 
or Mab ~ —21.0, and assuming that the other members 
of the same group are identical in nature to HUDF-JD2 
except for being fainter, the total power turns out to be 
between 4.2 and 6.7 times higher than HUDF-JD2 itself, 
depending on whether the LF extends just to L*/10 or 
down to L*/1000, respectively. As a consequence, the 
global effect for reionization would correspond to that of 
a source density ^ 5.4 ± 1.2 higher, i.e. 0.86 arcmin"^ 
(open circle in Figure 2). If this is the case, then it follows 
that HUDF-JD2, together with its lower luminosity com- 
panions, could easily produce enough Lyman-continuum 
radiation to reionize the IGM in their neighborhood un- 
der most reasonable circumstances {e.g., escape fraction 
as low as 0.2 and clumping factor as high as 20). Cor- 
respondingly the electron scattering optical depth for our 
reference model, in which reionization starts at zi — 13 
and is complete at Z2 = 10.7, turns out to be Tg = 0.13. 

3. DISCUSSION AND CONCLUSIONS 

We have considered the possible effects the source 
IIUDF-JD2 on the reionization of the Universe. We find 
that IIUDF-JD2 may indeed be capable of rcionizing its 
portion of IGM (i.e. the gas within a volume defined by 
the subtended solid angle of the HUDF-NICMOS field and 
the redshift of complete reionization, Z2 ^ 9 — 15), possibly 
with the help of unseen fainter companions, starting the 
process at a redshift as high as zi 2± 15±5. 

Regardless of whether HUDF- JD2 is able to account for 
the entire reionization of its cell of Universe, it is certain 
that it had the power to appreciably ionize it at high red- 
shifts. Even with no reddening correction and not allow- 
ing for faint companions, it appears that HUDF-JD2 alone 
should be able to ionize either 20% of its cell if the mat- 
ter distribution is homogeneous, or ionize it at an aver- 
age 20% level (i.e. the gas in the cell is 80% neutral) if 
there are density inhomogeneities in different directions . 
Mobasher et al. (2005) estimate that active star formation 
in HUDF-JD2 lasted less than 100 Myr, which implies a 
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redshift interval of Az 2 if the galaxy formed at z ^ 12 
and Az 6 if the galaxy formed at z ~ 20, during which 
HUDF-JD2 was ionizing its cell of the Universe at a ~ 20% 
level, at least. More generally, if HUDF-JD2 represents 
the "typical" source present in any 2.5' x 2.5' field on the 
sky, then this type of sources could readily ionize the IGM 
at levels and redshifts as high as required to account for 
WMAP findings (re = 0.17 ±0.06 at z= 17 ±5; Spergel et 
al. 2003). 

It is interesting to note that the four i'-dropout sources 
detected and studied by Eyles et al. (2005) in the 
HST/ACS GOODS images of the Chandra Deep Field 
South, may indeed be representatives of the "second 
brightest" galaxies in the ladder of the luminosity function. 
Actually, their measured magnitudes are about 1.5-2 mag 
fainter than measured for HUD-JD2. All of these sources 
are found to have z~ 5.8, stellar masses 2 — 4 x lO^^M©, 
ages of 250 — 650 Myr and implied formation redshifts 
Zf - 7.5 - 13.5 (Eyles et al. 2005). The presence of 
HUD-JD2 together with these additional bright sources 
in the Chandra Deep Field South (by design one of them, 
namely SMB03#1, falls within the HUDF area) appears to 
lend support to the possibility that the reionization of the 
Universe be dominated by massive galaxies that, rather 
unexpectedly, were able to form during the early evolu- 
tion of the Universe. 

Although the presence of a cluster of galaxies and the 
assumptions about their properties are somewhat specu- 
lative, we would like to stress that these are predictions 
that will be possible to confirm or exclude once the James 
Webb Space Telescope (JWST; e.g., Stiavelli et al. 2004), 



currently due to be launched in mid-2013, will come into 
operation. According to the model calculations by Pana- 
gia et al. (2003, and in preparation), sources like HUDF- 
JD2 and its analogs at higher redshifts will be easy to 
study spectroscopically because in their active star for- 
mation phases, they will have intrinsically high luminosi- 
ties providing fluxes as high as 150 nJy in the near and 
mid infrared. Fluxes at these levels are expected to be 
detectable with the JWST-NIRSpec spectograph at res- 
olution of ^ 1000 with a S/N ratio of 10 in 10''' seconds 
exposures. With the advent of JWST it will be possible to 
fully characterize these early Universe sources, thus clari- 
fying once and for all the process of galaxy formation and 
the reionization of the Universe. 
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Fig. 1. — Effective ionization parameter ^1^^ =< N^yc > / (yLi,)mm as a function of starburst age for different mctallicities (indicated 
as metal abundances by mass). For comparison, the values of the ionization parameters for black-bodies with temperatures of 50,000 and 
100,000 K are also shown. The unshaded area corresponds to the interval within which the age of HUDF-JD2 is likely to lie. 
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Fig. 2. — Surface density of possible reionization sources as a function of the apparent AB magnitude for identical sources with ionization 
parameter e|J^ = 9.6 x 10^^ erg~^. The lower solid line represents the minimum surfax;e brightness model, in which the ionizing radiation 
escape fraction / and the clumping factor C axe both unity, while the upper solid line represents the global metallicity constraint Z < O.OIZq 
at z = 6.5 (see SFP). The dashed lines represent models with {f,C) = (0.2, 1), (0.2,20), and (0.1,30), respectively, characterized by a lower 
ionizing efficiency. The position of HUDF-JD2 on the basis of its observed fluxes is denoted with a dot, whereas the triangle and the circle 
denote the positions after correction for dust extinction and after allowance for unseen companions, respectively (see text). 



